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Abstract 

The motion of a muon in two centers coulomb field is one of the interesting problems 
of quantum mechanics. The adiabatic expansion method is powerful approach to study the 
muonic three-body system. In this investigation the three-body problem is studied for short- 
range interactions. Bound states and energy levels of this system were calculated and compared 
with their Born-Oppenheimer method counterparts. The obtained results are in good agree- 
ment with the previous calculations. 

M.S.C. 2000: 36.10.-k, 21.45.+V 

Key words: three-body system, adiabatic expansion method, bound states, energy levels, 
Born-Oppenheimer method. 

1 Introduction 

A negative muon is a Icpton of the second generation with mass number about times heavier than 
that of electrons, and has a finite lifetime of = 2.197fisec. Nuclear fusion reactions can be 
catalyzed in a suitable fusion fuel by muons [1-3]. Energetic negative Muons, after stopping in 
a hydrogen isotope H2/D2/T2 mixture, fuse with different hydrogen isotopes to form the muonic 
atoms pfi, dfi and tfi in excited states. The size of a muonic atom is about 209 times smaller than 
a normal atom. After a sequence of cascade transitions lasting about 10~^^sec at liquid hydrogen 
density [L.H.D « 4.25 x atoms / cm?) , muonic atoms are formed in the ground IS'-state with 
kinetic energy in the range of 10^'^ — IQ^eV , depending on the temperature of the target and the 
prehistory of cascade transitions [4,5]. Muonic molecular formation may take place in collisions of 
the muonic atom in its ground state with ordinary hydrogen molecules. 

The crucial process is the resonant formation of the muonic molecule dtfj, into the loosely bound 
rotational-vibrational state J = 1, v ^ \ (lO^'* — lQ~^^sec). This state quickly de-excites to J = 
levels, from which the two nuclei fuse 10^^^). Usually the muon is free to start the next cycle, 
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but occasionally (with probability uJs — 0.555% in a D-T target at density (j) = 1.2 L.H.D [6]) it is 
removed from the active cycle by "sticking" to the helium nucleus. Resonant formation of the ddfi 
molecule at very low temperatures was observed in solid and liquid D2 targets [7,8]. Energy levels 
of muonic molecules in muon catalyzed fusion (fiCF) systems have been extensively studied with 
various calculational methods [9-12]. The solution of the non-relativistic Schrodinger equation 
for the three-body system was found using two different methods: by: the adiabatic expansion 
method and Born-Oppcnhcimer method. The Born-Oppcnheimer approach assumes the nuclei 
to be infinitely heavy with respect to the negatively charged particle. It should be kept in mind 
that the following Born-Oppenheimer approach is the simplest solution to the three-body coulomb 
system. This approach is expected to be a poor approximation for calculations of muonic molecule 
eigenvalues. In the present paper, we calculate binding energies of the bound states of the ddfi 
muonic molecule using the adiabatic expansion method. 



2 Theoretical calculations 

The exact Hamiltonian that describes muonic three-body system can be shown by following rela- 
tion: 

1-9 "^2 ^2 -^l-^fi ^2^j.L -^1-^2 

^ = '2^^^'^' " 2^^''- " 2;^^^" " Jr,, - Ril " Jr^ - R2I ^ JRi - R2I 

where 1 and 2 denote the two nuclei, their position is given by Ri and i?2, and the muon coordinate 
is r^. The center of mass coordinate Rcm is given by 

TOiRi + m2R2 + m^Vi_, . 
t^CM - ; ; 



It is convenient to define Jacobi coordinate r and R as follow: 

Ri + R2 
r = fi, 



(3) 



where R is the intcrnuclcar coordinate and r is the muon coordinate to midpoint between the 
two nuclei. In these coordinates {R,r), the Hamiltonian denoted by equation (1) is change to the 
following operator: 

where 

Mt = mi +1712 + rui^i (5) 
11 1 
M mi m2 
1 1 1 



m mi + 7712 mp 

mi - m2 

A — 

mi + m2 
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After separation of variables, the non-relativistic Hamiltonian in units of e = ?i = = 1, can be 
given by 



where 

^^i=-^V?-fl-f^ (7) 

2771 T"! f2 

where 6 represent the five dimensional variable. We use the set 6 = (0, ^, 77, (/?) where (0, 4>, (/?) 
define the Euler rotation specifying the body-fixed frame with its unit vectors to coincide with the 
principal axes of the inertia tensor of a three-body system. The hyperspheroidal coordinates ^ and 
r] are easily expressed by the muon- nucleus distances ri, 7'2 and the intcrnuclear distance R, 

(8) 

The three-body Hamiltonian (6) commutes with the total angular momentum operator for the 
three particle system, J, its projection on z-axis, J^, and the total parity operator, P(R — > 
-R, r — > -r) . Eigcnfunctions of the Hamiltonian in the total angular momentum representation 
reads: 

vI/%,(R,r)= ^ Fi^{R,tv)DiU<^',Q,^) (9) 

m— — J 

Adiabatic expansion of radial function, ^, //) is usually written in the form: 

00 N — l 00 
Af=l 1=0 1=0 •'^ 

where Xi^{R) describe relative motion of the nuclei. Let us consider the Wigner function, D'^fj^i^, &, (fi) 
which is the eigenstates of J^, Jz and R.J/R with the eigenvlues J(J -t- 1), M and 771 [13]. It can 
be transformed under the inversion as follow: 

PDi,^{<^>, e, ^) - Di^{^ + 7T,n-e,n-^) = (-l)'^-"Z?)f,^_„($, 9, ^) (11) 

If 771 7^ 0, the resultant Wigner functions would be different, and the angular functions consist both 
even and odd combinations. It is convenient to specify these combinations as follows: 



P = ±(—1)'^ is the eigenvalue of the parity operator: 

PDifn^-pDil^ (13) 

The functions presented in equation (12)(in bracket) are satisfy the following orthonormality con- 
dition: 



p2n 

iinQdQ / dip 



KL. 0, A K'..' 0' ^) = ^Jf V ^MM' ^nm' (14) 
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If m = 0, both the Wigner functions in (11) are reduced to the ordinary spherical function 
Yjm{Q, ^) so that the dependence of (p disappears and the angular functions satisfying the condi- 
tions (13) and (14) are: 

rip ,^ o ,_,_Yjm{Q,'^) 



D 



($,e,^) 



(15) 



In this case the parity is unambiguously specified by the quantum number J: p = +( — 1)'^. So, 
our basis functions have the following structure: 



R 



(16) 



The wave functions ^'][]'^^^(i?, 0, $, ^, 77, (p) describing reactions /i/i + /i, h = (p,d,t) can be de- 
composed over the solutions ipjm[^^ri; R) of the Coulomb two-center problem. ipjjni£,,r]j R) is 
complete set of solutions of the Coulomb two-center problem, therefore 



HitP,{^, 77; R)F{^) = E,{R)ij,{^, 77; R)F{^) 



(17) 



describing the muon motion around fixed nuclei separated by a distance R. Ei (i?) is the energy of 
a muon in the state i as a function of R. Here we show how to separate the variables through the 
use of the ellipsoidal (or, prolate spheroidal) coordinates 



1)(1 — lif)cOSLp 



(18) 



y = - 1)(1 -■rf)sirup 



Note that the coordinates ^, rj and ip arc orthogonal, and we have the first fundamental form 



ds^ = dx^ + dy^ + dz^ = hpf + h^drf + h^dip^ 



(19) 



where 



dx 
dx 

dx 

dip 



dy 
drj 

dy 

dip 



dz 
dl 
dz 
di] 

dz 
dp 



_ R^ 
4 

- ^ 
4 



e 
e 



(20) 



Thus 



1 



d 



d_ 

9C 



hr,h^ d 



d_ 

drj 



d / h^h^ d 
drj \ hri drj 
d 



d / h^hn d 
dip \ dip 

e - if 9' 
(e2-l)(l-r;2)5^2 



(21) 
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Note that through the coordmate transformation (18), we have 



r2 



(22) 



Writing the wavcfunction as R)F{'p) ~ G{£,)H{7])F{ip) and changing the variable to spheroidal 

coordinates, equation (17) can be separated into following three one-dimensional equations: 



d_ 
dr} 



d( 



drj 



-A 



-A 



-Tn^Fiif) = 



(iqrj + q^rf 



e-i 

^2 



1 



G(0 = 
H{i]) = 



where 



a 



R , , 

— [Zi + Z2) 

q 

R 



(23) 

(24) 
(25) 

(26) 



^2) 



Note that A and q are unknown parameters and should be obtained from (24) and (25) as eigenval- 
ues of the coupled system. Once A and q are obtained, then E can be obtained from q^ ~ —B?E/2. 
By substitution of expression (16) into the Schrodinger equation with Hamiltonian (6) and after 
averaging over spherical angles (0, $) and the muon state, one obtains the radial equation 



1 d^ 



e~VB-o{R)~U^{R) 



J{J+l) 



xl{R) = 



(27) 



where e = E—Ei{oo) is the collision energy and E is the total energy of the system and Ei{oo) is the 
ground state energy of muonic atom. VB-oiR) — E — E{oo) + ^^ is the potential corresponding to 
the Born-Oppenheimer (B-0) approximation and Ui{R) = (i| (iJ — Hi — \i) is the adiabatic 
correction. The adiabatic potential VAd{R) is: 



VAdiR)^VB-o{R) + U,iR) 



(28) 
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Table 1: Binding energies {eV) of the states (J, j^) for the (dd/i) muonie molecule. 





Ad (This work) 


Ref. [14] 


Ref. [15] 


(0,0) 


325.06 


325.0735 


325.070 540 


(0,1) 


35.79 


35.8436 


35.844 227 


(1,0) 


226.62 


226.6815 


226.679 792 


(1,1) 


1.73 


1.97475 


1.974 985 


(2,0) 


86.20 


86.4936 
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Figure 1: Adiabatic potential curves, VAd(J = 0, 1 and 2), corresponding to d — d — fi system. 

The adiabatic potential VAd{R) for the {ddji) muonie molecule is calculated in the adiabatic expan- 
sion method. The adiabatic potential curves and qualitatively similar for each of muonie molecules 
and are displayed for {ddjj) in Figure 1. Results of the calculations of binding energies of the bound 
states (J, v) of the (ddfi) muonie molecule are compared with the results of the other methods used 
in Refs. [14,15] and are given in Table 1. 



3 conclusion 

In this paper, we have presented the calculation of the binding energies of the symmetric muonie 
molecular ion. Early calculations were unable to demonstrate the existence of the crucial weakly 
bound states. For example, in the Born-Oppenheimer (fixed nuclei) approximation the state is 
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much too bound, but if adiabatic corrections are included it is not bound at all. The results 
are obtained by the three body Hamiltonian in the adiabatic expansion method. The calculated 
binding energies are in good agreement with the previous calculations by other authors using 
different methods. 
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